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The specific heat and susceptibility data for 3 He on Grafoil are analyzed in the 
framework of the Landau Fermi liquid theory. The dominant interaction between 
3 He quasiparticles is found to be in the p-wave channel for most experimentally 
accessible areal densities of 3 He. This interaction is attractive and gives rise to the 
p-wave transition temperature which for moderate areal densities is estimated to be 
on the scale of several millikelvin. The relevance of these results to the anomaly in 
the specific heat observed at Tfc = 3.2 mK is discussed. 
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Experiments on monolayer films of 3 He absorbed on Grafoil [|T],0,[3[ and on 3 He- 4 He mix- 
ture films on Nucleopore substrate 0|| provide an opportunity to study in detail the prop- 
erties of 2D Fermi liquids. Interest in this subject has grown since the validity of the 
conventional Fermi liquid theory in 2D was questioned in the context of the high-T c super- 
conductivity ||. In the present approach, we focus on the 2D Fermi liquid at sufficiently 
low densities where perturbative calculations ]7| do not show any divergencies which might 
signal the breakdown of the Fermi liquid description. Accordingly, we assume that the Fermi 
liquid picture is valid for 3 He in two dimensions. 

The most remarkable property of bulk 3 He is a superfluidity with nonzero angular mo- 
mentum I = 1 ||. However, in "surface" 3 He the superfluid transition has not yet been 
observed. In what follows, we analyze the possibility of superfluidity for 3 He on Grafoil on 
the basis of information which can be inferred from experiments on specific heat |l[ and 
magnetic susceptibility ||. We find that for most experimentally accessible densities, this 
system cannot be described by a momentum independent (i.e., s-wave) interaction. Instead, 
the dominant interaction component is in the p-wave channel. The corresponding scattering 
amplitude is attractive so that for moderate densities one might expect to get relatively 
high T c ~ lOOmK. However, our calculation yields an anomalously small prefactor in the 
expression for T c in 2D which shifts the transition down to the millikelvin region. The cal- 
culated value of T c is reasonably close to 3.2 mK, where the specific heat anomaly in 3 He on 
Grafoil has been observed We argue that this anomaly may correspond to the onset of 
superfluidity. 

We start with a brief review of the properties of dilute Fermi liquids 0. At small 
densities, the s-wave component of the scattering amplitude is dominant. The perturbative 
expansion holds in powers of the product of the scattering amplitude and the density of 
states at the Fermi surface, iVo- In three dimensions N oc pp while the s-wave scattering 
amplitude (scattering length) is normally of the same order as the range of the interaction 
potential, r . In two dimensions, the density of states N is independent of p, but a low 
density expansion is still possible because the scattering amplitude in 2D tends to zero 
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logarithmically as p — > |1(J. Accordingly, the expansion parameter for the 2D problem is 
g(p) = l/log(p /p), where p ~ r o 2 - When p <C p , the expressions for the effective mass 



m* and magnetic susceptibility x m 2D are fllTf : 



— = 1 + ^1 + 2^, 

m 3 

y 1 + Ff m* / 9 , \ -i 

X ' ~ — (l-2^ + Vln2) , (1) 



Xz 1 + F a ™3 

where and %3 refer to an ideal gas of 3 He atoms. 

The experimental data on the density dependence of the specific heat and magnetic 
susceptibility are presently available for 3 He on Grafoil |l],@,[§ and for 3 He- 4 He films on the 
Nucleopore substrate ||,||. The latter case is more difficult to analyse because a 3 He atom 
occupies a surface bound state on top of 4 He and its hydrodynamic mass is substantially 



larger than the atomic mass of 3 He due to the interaction with underlying 4 He layers [12 
In what follows, we concentrate solely on the properties of the 3 He film on Grafoil. The 
experimental results for the effective mass |3J along with theoretical predictions for the s- 
wave and p-wave scattering amplitudes are presented in Fig.[l|. Although the higher order 
terms in g may be important for larger densities, it is clear from Fig.|TJ that the (logarithmic) 
density dependence of g is too weak to account for the fast increase of the effective mass 
as density increases. The above discrepancy signals that in the experimentally accessible 
region of densities the I = harmonic does not overshadow the higher angular momentum 
components. In general, in this situation all harmonics should have equal strength. However, 
we found that the plot of m*/m versus p is well described by a simple fit which involves only 
the p-wave component of the scattering amplitude (Fig.[l|): 

Ti-nd P V 1 (2) 

m V 0.062 atoms/A 2 / 

The nearly linear behavior of m/m* as a function of p which is seen in the experiments |l| 
even at relatively small p indicates that for 3 He absorbed on Grafoil the p-wave component 
of the scattering amplitude is anomalously large so that it overshadows the contribution 
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from the s-wave channel in nearly the entire experimentally accessible range of densities. 
This fact inspires us to reexamine the low-energy expansion for 3 He on Grafoil. 

To proceed with the theoretical description, consider first the case when the Born ap- 
proximation is valid, i.e. the scattering amplitude coincides with the spin independent 
interaction potential U(k' — k). Simple algebra then yields F a = — (m*/27ih 2 )U(k' — k), 
F s = {m* /27rh 2 )(2U(0) — U(k' — k)). Here both k' and k are on the Fermi surface, so that 
F a,s depend only on the angle 9 between them. The interaction potential can be expanded 
at low densities as U(k' - k) = U(0) + A(k' - k) 2 + ... = U(0) + 2Xp 2 F (l - cos#) + 0{p%). 
This immediately yields the following expressions for the effective mass (same as @) and 
spin susceptibility: 



^ ~ l-Xm 3 p, ~ ^ (l-2g-2\m*p) . (3) 

Leaving the detailed comparison with the experiment to the study of the effects beyond 
the Born approximation, we merely conclude at this point that in order to account for the 
increase of the effective mass, A should be positive so that the pairing interaction in the 
p-wave channel is attractive |]T3. The coupling constant for the p-wave pairing is of the 
order of 1 — m 3 / m* . Since the experiment shows that the effective mass may well exceed the 
bare mass and ep is nearly IK, it is not clear a priori why no superfluid transition has been 
observed down to the millikelvin temperature region. To address this issue, we now calculate 
explicitly the transition temperature in a 2D Fermi gas with p-wave attractive potential and 
show that the prefactor in T c is anomalously small in two dimensions. 

A way to calculate the prefactor in T c in a weak coupling approximation is to start with 
perturbation theory and collect all second order contributions to the pairing vertex which 



come from the integration over momenta far from the Fermi surface |14|]. The renormalized 
vertex should then be substituted into the Cooper channel and the integration within the 
ladder should be restricted to a region close to the Fermi surface where the logarithmic term 
is dominant. By carrying out the above procedure, one obtains an equation for T c [ H ]: 

l = g(l + ag)ln e 7 f + 0(g 2 ), (4) 
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where g is a coupling constant and a is a numerical factor. Solving (H), one gets T c = 
€p exp(— 1/g), where ep = e^exp(a). 

Although for any interaction strength the weak coupling approximation is valid 
at sufficiently low density, the renormalized p-wave vertex in the Cooper channel 
T(p, — p; p', — p') = T(9) (cos 9 = pp'/Pf) generally cannot be expressed in terms of a 
single parameter because the total p-wave scattering amplitude is an unknown nonuniversal 
quadratic function of the momenta. In other words, unlike for the s-wave case, here one 
can not substitute the same scattering amplitude into all vertex functions in the second 
order diagrams which contribute to T(6). In view of this, we first perform the calculation 
assuming that the Born approximation is valid. Then we reconsider the problem by taking 
into account some of the effects beyond the Born approximation. 

There are four second order diagrams which contribute to T(6). Three of them are from 
the zero sound channel (Fig.0a-c) while the fourth is from the Cooper channel (Fig.|2]d). The 
evaluation of the zero sound diagrams is lengthy but straightforward. We take advantage 
of the fact that in the region of experimental interest the s-wave component of the interac- 
tion is considerably smaller than the p-wave component and neglect U(0), i.e. substitute 
U(kiak 2 f3; k 3 7k 4 cr) = A(k x — k 3 ) 2 5^5^ into the vertices of Fig.|2|. The calculation then 
yields 

97m\ 2 p F 

i a-c = lr t2 CQS A \°arfOp<r + <W<W ■ (5) 

i&irn 

The last diagram (Fig.pjd) contains both log(ep/T c ) and the contribution to the prefactor, 
and also contributes to the vacuum renormalization which transforms the interaction po- 
tential into the scattering amplitude 0. The vacuum renormalization has to be subtracted 
from the diagram of Fig.|2|d. This procedure eliminates ultraviolet divergence in the theory 
leaving ep as the only dimensional quantity: 
r _ flmh 2 ^\2e Fl 
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cos 6 (5 ai 8p a + 8 a<7 8p 7 ) , (6) 

where In 7 is the Euler constant C ~ 0.577 and fi is the p-wave scattering amplitude, which 
in the Born approximation is equal to —2Xp'j 7 /h 2 . 



Expressing ([5]) in terms of f\ and combining it with (|6]), we get the renormalized pairing 
interaction in the p-wave channel T = T a _ c + which includes both leading (logarithmic) 
and next to leading order terms in f%. The instability criterion then yields (fi < 0): 



While 27/7? m 1.13, the factor exp(— 112/15) ~ 5.72 • 10~ 4 reduces T c in (|7j) by more than 
three orders of magnitude. As a result, even if the coupling constant g p = m|/ 1 |/47r ~ 1, 
the critical temperature T c ~ 10~ a €f- 

The results above were obtained in the Born approximation, i.e. under the assumption 
that the Born parameter u ~ Xm^/AnrQ <C 1. The validity of the Born approximation 
is, however, questionable for 3 He on Grafoil. Indeed, it follows from that in the Born 
approximation the linear in density term in Fq is two times larger than that in F{. However, 
the experimental results [1],[J give a much smaller value for this ratio. This suggests that 
the Born parameter is, in fact, not small and hence the density-independent ("vacuum") 
corrections to the Fermi liquid parameters are important. We calculated the leading vacuum 
p-channel corrections in the symmetrical gauge and obtained 



Since A should be positive in order to account for the increase of m*/m 3 with the density, 
vacuum corrections increase the linear term in Ff and reduce that in Fq which narrows 
the gap between theory and experiment. A fit to the experiment gives u ~ 1—2, that is 
the vacuum corrections are, indeed, strong (the uncertainty in u is related to the s-wave 
contribution to Fq which is difficult to estimate precisely). 

The next step would be to calculate T c beyond the Born approximation. However, 
we already mentioned that in order to solve the problem one needs to know what the 
renormalized scattering amplitudes are for all combinations of the momenta relevant to the 
diagrams of Fig.|^. Such calculation is lengthy and not very informative because it is not clear 






(8) 
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whether one can deal with only the leading term in u in the region of parameters relevant to 
the experiments. Because of this complication, below we use a more qualitative approach and 
just take into account the fact that the spin structure of the p-wave interaction potential does 
not survive the effect of vacuum renormalization, i.e. the total p-wave scattering amplitude 
has both spin-independent and spin-dependent parts even if the initial interaction was spin- 
independent. Accordingly, we model the effect of vacuum renormalization by introducing 
an effective potential which satisfies the Born approximation condition and has both spin- 



independent and spin-dependent parts [fL5 



U ef f = U (k' - k)^^ + U(k' - k)a ai a, 3a . (9) 

The low density expansion of U is U(9) = {7(0) + 2\p 2 F (l — cos 6) + 0(p F ). This effective 
potential reproduces measured m*(p) and x(p) quite well for £ = A/(A + A) ~ 0.3 — 0.4 (note 
much smaller uncertainty in ( than in u) . 

The calculation of T c with U e ff proceeds along the same lines as above. We skip the 
intermediate calculations and present only the result: 

2 7 e F / 112 56 308 2 \ f l\ 

The p-wave coupling constant now is 

9 P =^(X + X). (11) 

27111 

As it turns out, ( = 0.3 — 0.4, inferred from the fit to the experiment, corresponds to the 
prefactor ep = 0.03 — 0.04 ef- This value is larger than our earlier estimate though is still 



considerably smaller than the p-wave result in three dimensions, O.lep ||16|| . It is worth 
mentioning here that for bulk 3 He the value of the prefactor inferred from the measured T c 



was found to be rather insensitive to the particular form of the interaction |I7 
In order to estimate T c , we express the coupling constant g p as 



For small densities, m* m and T c is exponentially small. However, for moderate densities, 
m* is significantly larger than the bare mass and the coupling constant saturates at the value 
(1 + 2C) _1 « 0.5—0.6. Substituting this value of g p into (|T0|), we get T c ~ 5 — 7mK which is 
of the same order as T k pa 3.2 mK where the anomaly in the specific heat has been observed. 
Although the quantitative agreement with the experiment is unanticipated because of the 
approximate nature of the theoretical considerations, our results indicate that the dense 
"surface" 3 He on Grafoil may become superfluid in the experimentally accessible temperature 
range. Note that there is no true off-diagonal long range order in two dimensions, but the 
superfluid density is finite below the transition |18|,|19 |. The actual Kosterlitz-Thouless 



transition temperature does not differ substantially from the "mean-field" T c which we have 
calculated 

It is argued in that the specific heat anomaly in 3 He on Grafoil at Tj, is an intrinsic 
property of the fluid monolayer. The calculated value of T c flip]) is comparable with T k , 
but, unlike T k , it significantly drops down at low densities. We suggest that the density 
independence of T k observed in the experiments may in fact be a result of a phase separation 
which accompanies the superfluid transition, that is for arbitrary density of 3 He there is an 
energetically stable superfluid dense phase of 3 He below T k . This tentative scenario explains 
both the lack of density dependence in T k and a rapid decrease of the specific heat below 
T k . 

To summarize, in this paper we have analysed the specific heat and susceptibility data 
for 3 He on Grafoil. We found that the p-wave component of the interaction between 3 He 
quasiparticles is dominant for all experimentally accessible areal densities of 3 He. This 
interaction is attractive and gives rise to a p-wave transition temperature which for the 
dense "surface" 3 He is estimated to be in the millikelvin temperature region. We suggest 
that 3 He on Grafoil may be superfluid below 3.2 mK where the anomaly in the specific heat 
has been observed. 
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FIGURES 

FIG. 1. The density dependence of the effective mass for 3 He on Grafoil. The experimental 
data Q is fitted using the interaction potential with dominant s-wave or p-wave amplitudes (dotted 
and dashed lines, respectively). 

FIG. 2. Second order diagrams which contribute to the pairing vertex. The wavy line represents 
the interaction potential. 
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